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SUMMARY 

In  a  recent  work,  we  presented  a  new  technique  for 
analyzing  the  asymptotic  nature  of  solutions  of  linear 
differential  equations  and  linear  differential-difference 
equations.  We  mentioned  there  that  the  same  method  is 
applicable  whenever  a  certain  integral  representation, 
essentially  a  convolution  transform,  exists  for  the  solution. 
The  purpose  of  this  note  is  to  initiate  the  application  of 
the  method  to  the  study  of  partial  differential  equations 
possessing  thlB  property. 

Here  we  consider  the  parabolic  equation 

t 

ut  "  uxx  +  >(x)a(t)uf 

with  the  boundary  conditions  u(0,t)  *  u(l,t)  »  0,  t  >  0. 

In  subsequent  papers,  we  shall  apply  the  method  t>o  other 
types  of  functional  equations.  For  the  sake  of  simplicity, 
assume  that  a(t)  -  (t  +  l)-1  or  (t  +  l)-2,  and  that 
j^(x)  is  a  continuous  function  over  0  <  x  <  1. 
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ASYMPTOTIC  BEHAVIOR  OP  SOLUTIONS 
OF  LINEAR  PARABOLIC  EQUATIONS 

Richard  Bellman,  The  RAND  Corporation 
Kenneth  L.  Cooke,  Pomona  College 

1 .  INTRODUCTION 

In  a  recent  work\  we  presented  a  new  technique  for  analyzing 
the  asymptotic  nature  of  solutions  of  linear  differential 
equations  and  linear  differential— difference  equations.  We 
mentioned  there  that  the  same  method  is  applicable  whenever  a 
certain  integral  representation,  essentially  a  convolution  trans¬ 
form,  exists  for  the  solution.  The  purpose  of  this  note  is  to 
initiate  the  application  of  the  method  to  the  study  of  partial 
differential  equations  possessing  this  property. 

Here  we  consider  th6  parabolic  equation 

ut  -  uxx  ♦  ^(x)a(t)u,  (1) 

with  the  boundary  conditions  u(0,t)  -  u(l,t)  ■  0,  t  >  0.  In 
subsequent  papers,  we  shall  apply  the  method  to  other  types  of 
functional  equations.  For  the  sake  of  simplicity,  assume  that 
a(t)  »  (t  +  l)-1  or  (t  +  l)-2,  and  that  /(x)  is  a  continuous 
function  over  0  <  x  <  1.  Then  we  can  establish 
Theorem.  For  each  positive  Integer  n,  let 

cn  -Z1  ^(x1)sln2nvx1dx1 .  (2) 

If  a(t)  -  (t  +  l)-2,  for  each  n  there  le  a  solution  of  (1), 
satisfying  the  stated  boundary  conditions,  with  the  asymptotic  form 

2  2  2  2 

u(x,t)  »  e-^'  T  tsin  mrx  +  Ofe-0  v  ^t  *)  (.3) 

as  t  — ►  oo  . 

If  a(t)  «  (t  +  l)_i,  there  is  a  solution  of  the  form 
u(x,t)  -  e-11  *  tsin  nvx  +  ofe-11  *  t) 
if  cn  <  0,  and  of  the  form 


(*) 
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22.  2c  22.  2c 

u(x,t)  -  e-*1  H  nsin  mrx  +  H  )  (5 

If  cn  >  0.  The  Implicit  constants  In  the  o— relations  are  uniform 
In  x  for  0  <  x  <  1 . 

It  will  be  clear  from  what  follows  how  more  general  forms  of 
a(t)  are  handled. 

i 

2.  CONVOLUTION-TYPE  INTEORAL  EQUATION 


It  is  easy  to  see  that  the  solution  of 


ut  •  u*x  + 

u (x,0)  -  0,  0  <  x  <  1, 

u(0,t)  -  u(l,t)  -  0,  t  >  0, 


may  be  written  In  the  form 

..  /)  t 


U  m,Jo  efc  k^x,xl#t  “  t1)f (x1,t1)dx1dt1 


where 


QD  2  2 

k(x,x.,8)  *  2  ^  e  Baln  mrx  sin  mrx^ 
n-i  2 


-  e^(*(x  -  x1)/2,e~J*  t)  -  e^(ir(x  +  x1)/2,e  T 


In  the  notation  of  Whittaker  and  Watson  . 

Observe  the  convolutary  nature  of  the  Integral . 


3.  NON-ABSOLUTELY  INTBORABLE  a(t) 

Since  the  case  where  a(t)  Is  absolutely  lntegrable  over 
[0,co]  Is  rather  simple  to  treat,  let  us  concentrate  on  the  more 
Interesting  and  difficult  case  where  this  is  not  so.  Let  us  take 
a(t)  »  (t  +  l)-1  as  an  example,  n  ■  1,  and  c^,  as  defined  by 
(1.2),  to  be  nonnegative. 

Applying  the  result  of  (2.2)  to  (1.1),  with  f(x,t)  ■  /(x)a(t)u, 
we  obtain  a  linear  integral  equation  which  we  write  In  the  dis¬ 
jointed  form 
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u(x,t)  -  e  T  sin  vx  +  p(x,t)  (1) 

i.  .  •«  — 1 1  (t— t,  ) 

♦  2  sin  vx  J'  J*  e  sin  wx^fx^ )  (t.^4-1 )  u(x^ ,  Jdx^dt.^ 


where 


p(x,t)  mc/>t(/)l*1  (x#x1,t-t1)/(x1)  (t^l)  1u(x1,t1)dx1dt1, 

oo  ^2  2 

k-L(x,x1,t)  «  2  ^  e”"*1  *  sin  nvx  sin  mrx^ . 

n-2 

TYie  expression  In  (1)  suggests  that  we  introduce  a  new 
variable 

v(t )  -jO1  ^(x1)sin  irx^u  (x1,t1)dx1, 
and  write 

2 

u(xft)  ■  e_T  tsln  vx  +  p(x,t)  c 

t  — ir2(t-t,) 

+  2  sin  vXj/1*'  e  1  (t1  ♦  1)  ^(t^Jdt^. 

Multiplying  by  /(x)sln  nx  and  integrating  over  x,  we  derive 
the  relation 

2.  .  -v2 ( t~t,  )  '  , 

v ( t )  -  Cxe  1  4-  2ci^/?t  0  1  (tx  +  1)  v(t1)dt1  4-  q(t), 

where 

q(t)  •  p(x,t)^(x)sln  txdx. 

Following  the  procedure  used  in  \  we  replace  (5)  by  a  more 
tractable  Integral  equation.  Differentiating  (5),  we  obtain 

v»  -  q*  -  (-  v2  +  2c1/(t  4-  1) )  (v  -  q)  +-  20^(1  4-  1), 


,-l 


(2) 


(3) 


(*) 


(5) 


(b) 


(7) 


which  yields  the  relation 

g  ^  ^  j 

V (t )  -  C.9,(t)  +  e  1  (q  +  1  )_1q(t1  )<lt1  +  q(t),  (6) 

where  s(t)  *0/,t  (—  +  2c^/(t^  +  l))dt^. 


P—1870 

4 


This  equation,  together  with  the  equations  relating  u,  v, 
p  and  q,  can  be  uaed  to  derive  the  stated  result.  The  same 
method  can  be  used  to  handle  more  general  functions  a(t). 
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